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Abstract: 

We analyze the implications for the status and prospects of supcrsymmctry of the Higgs discovery 
and the last XENON data. We focus mainly but not only on the CMSSM and NUHM models. 
Using a Bayesian approach we determine the distribution of probability in the parameter space 
of these scenarios. This shows that, most probably, they are now beyond the LHC reach . This 
negative chances increase further (at more than 95% c.l.) if one includes dark matter constraints 
in the analysis, in particular the last XENON100 data. However, the models would be probed 
completely by XENON1T. The mass of the LSP neutralino gets essentially fixed around 1 TeV. 
We do not incorporate ad hoc measures of the fine-tuning to penalize unnatural possibilities: such 
penalization arises automatically from the careful Bayesian analysis itself, and allows to scan the 
whole parameter space. In this way, we can explain and resolve the apparent discrepancies between 
the previous results in the literature. Although SUSY has become hard to detect at LHC, this does 
not necessarily mean that is very fine-tuned. We use Bayesian techniques to show the experimental 
Higgs mass is at ~ 2 a off the supersymmetric expectation. This is substantial but not dramatic. 
Although the CMSSM or the NUHM are unlikely to show up at the LHC, they are still interesting 
and plausible models after the Higgs observation; and, if they are true, the chances of discovering 
them in future dark matter experiments are quite high. 
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Physics, LHC Physics, Dark Matter Searches. 



Contents 



1. Introduction 1 

2. Set up 3 

2.1 The Bayesian approach 4 

2.2 Model-independent incorporation of naturalness criteria in the MSSM 5 

2.3 Priors in the initial parameters 7 

3. Physical observables and computational methods used in the analysis 9 

4. Results for the CMSSM 12 

4.1 Results with no dark matter constraints 13 

4.2 Incorporation of dark matter 17 

4.3 g-2 24 

5. NUHM 26 

6. Can still SUSY be natural? The impact of the Higgs discovery in the 
naturalness of SUSY 27 

7. Comparison with previous literature 33 

8. Conclusions 35 
A. Focus-Point behavior in the CMSSM and the NUHM 36 



1. Introduction 

The ATLAS and CMS collaborations have reported [1, 2] the strong evidence evidence of a 
Standard Model (SM) -like Higgs boson with rrih ~ 126 GeV. According to the most recent 
analysis presented in the Kyoto conference [3, 4], which combines the y/s = 7 TeV (L=5.0 
ftr 1 ) and y/s = 8 TeV (L=12.1 fb" 1 ) data, the "Higgs-boson" mass is 

m h = 126.0 ± 0.4 (stat) ± 0.4 (syst) GeV (ATLAS) 

m h = 125.7 ±0.4 (stat) ± 0.4 (syst) GeV (CMS). (1.1) 
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Adding the statistical and systematic uncertainties in quadrature, and combining further 
the ATLAS and CMS results, using the principle of maximum likelihood, one finds the 
following central value and uncertainty 



Beyond the crucial importance of this discovery by itself, this result has far-reaching 
consequences for well- motivated candidates for physics beyond the SM, such as Supersym- 
metry (SUSY), and in particular the Minimal Supersymmetric Standard Model (MSSM). 

Many works have explored these implications from different points of view, often re- 
lying on popular versions of the MSSM, such as the constrained MSSM (CMSSM) or the 
(slightly less restrictive) non-universal Higgs masses model (NUHM). In addition, some 
analyses incorporate as well the recent XENON100 constraints on fermionic dark matter 
[5, 6, 7, 8, 9, 10, 11]. Roughly speaking, the MSSM is in trouble, since a ~ 125 - 126 GeV 
Higgs generically requires rather large SUSY masses, which is in tension with the natu- 
ralness of the electroweak breaking at the correct scale. A big portion of the previously- 
acceptable CMSSM and NUHM parameter-space becomes excluded. Some of these ex- 
plorations present the results in the form of scatter plots, which illustrate the acceptable 
regions in slides of the parameter space, see e.g. [10]. While these analyses certainly reflect 
the situation, it is more difficult to tell from them the degree of tension of the CMSSM and 
the NUHM (or the general MSSM) with the experimental result rrih ~ 126 GeV or what 
is the relative probability of the various regions of the associated parameter-space. 

There are two main strategies to attack these questions: the frequentist and the 
Bayesian approaches. The frequentist one is based on the analysis of the likelihood func- 
tion in the parameter space: one performs a scan of the parameter space, evaluating the 
likelihood based on the \ 2 . This leads to zones of estimated probability (inside contours 
of constant x 2 ) around the best fit points. In doing so, the unplotted variables are opti- 
mized to obtain the best \ 2 . A nice feature of the frequentist approach is that it does not 
require to establish a prior to analize the parameter space. Essentially, it corresponds to 
a search of the points with best likelihood, and therefore gives relevant information about 
the goodness of different regions of the parameter space. However, a frequentist approach 
does not give necessarily the fair probability-map of the parameter space of a model. This 
is obtained using a Bayesian approach. 

As an example of this, let's take the notion of fine-tuning. Fine-tuning has to do with 
statistical weight and a frequentist analysis is based entirely on likelihood, i.e. the ability 
to reproduce the experiment, and thus cannot "see" the fine-tuning. E.g. a point (o a 
region) in the parameter space can present an optimal likelihood, but only after an extreme 
tuning of the unplotted parameters. This feature makes this point extremely unplausible 
(even if it is capable to reproduce the experimental data), but this does not affect its 
privileged condition in a frequentist approach. On the other hand, a correct Bayesian 
analysis incorporates automatically (without any ad hoc consideration) a penalization of 
the fine-tuned regions, since they have little statistical weight. 



m h = 125.85 ±0.4 
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Therefore the Bayesian approach offers a representative diagnosis of the health of a 
theoretical scenario and the future prospects. The main goal of this paper is to perform such 
Bayesian diagnosis in full detail, incorporating the impact of a possible Higgs at ~ 126 GeV, 
as well as the last XENON100 dark-matter constraints and other recent experimental data. 
We will also incorporate significant improvements over previous Bayesian explorations. 
They will allow us to avoid the arbitrariness associated to the choice of the ranges where 
the parameters live. We will show that this ambiguity is the origin of important apparent 
discrepancies between previous literature. From the analysis we obtain the most probable 
regions in the CMSSM and NUHM parameter space and study the capacity of the LHC 
and XENON of probing those regions in the future. 

Using Bayesian techniques, we will also explore in more general terms how the Higgs 
discovery affects the naturalness of SUSY, quantifying the tension between ~ 126 GeV 
and the supersymmetric predictions. 

In sect. 2 we establish our set up, reviewing the Bayesian approach and how it can 
be rigorously implemented in the context of the MSSM. We discuss also the choice of 
priors issue. In sect. 3 we list the observables used in the analysis in order to calculate the 
likelihood of the different points of the parameter space. We also expose the different codes 
and computational techniques used in the numerical calculations. Sect. 4 is devoted to the 
results for the CMSSM. We present them with and without dark matter constraints. We 
also distinguish between the possibility that the dark matter is purely supersymmetric or 
it has other components. We also quantify the discrepancy between present experimental 
data on g — 2 and the Higgs mass in the supersymmetric framework. This quantifiable 
tension is the reason why we have chosen not to include g — 2 data in the analysis until 
the theoretical situation becomes fully clarified. In sect. 5 we perform a similar analysis 
for the NUHM and present the results. In sect. 6 we study the issue of the naturalness 
of SUSY after the Higgs discovery, using Bayesian techniques. Sect. 7 is devoted to 
the comparison with previous literature on this subject and to discuss the reasons of the 
apparent discrepancies. The conclusions are presented in sect. 8. In addition we describe 
some meaningful theoretical aspects of the Focus-Point region in the CMSSM and the 
NUHM. 

2. Set up 

We will consider two versions of the minimal supersymmetric standard model (MSSM), 
namely the so-called constrained MSSM (CMSSM) and the non-universal Higgs masses 
model (NUHM). In the CMSSM the soft parameters are assumed universal at a high scale 
(Mx), where the supersymmetry (SUSY) breaking is transmitted to the observable sector, 
as happens e.g. in the gravity-mediated SUSY breaking scenario. The CMSSM parameter- 
space is then defined by the following parameters: 

{6i} = {m ,m 1/2 ,A,B,[i,s} . (2.1) 
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Here too, TO1/2 and A are the universal scalar mass, gaugino mass and trilinear scalar 
coupling; B is the bilinear scalar coupling; fi is the usual Higgs mass term in the superpo- 
tential; and s stands for the SM-like parameters of the model, i.e. essentially gauge and 
Yukawa couplings. All these initial parameters are understood to be defined at Mx- The 
NUHM is defined by the same parameters as the CMSSM but allowing for a different scalar 
mass for the two Higgs doublets, to#. 

2.1 The Bayesian approach 

The goal of the Bayesian approach is to generate a map of the relative probability of the 
different regions of the parameter space of the model under consideration (CMSSM and 
NUHM in our case), using all the available (theoretical and experimental) information. This 
is the so-called posterior or probability density function (pdf), p(0j|data), where "data" 
stands for all the experimental information and 9i represent the various parameters of the 
model. The posterior is given by the fundamental Bayesian relation 

p(0i|data) = p(data|0i) p(0;) — , (2.2) 

p(data) 

where p(data|#j) is the likelihood (sometimes denoted by £), i.e. the probability density of 
measuring the given data for the chosen point in the parameter space (this is the quantity 
used in frequentist approaches); p(9i) is the prior, i.e. the "theoretical" probability density 
that we assign a priori to the point in the parameter space; and finally, p(data) is a 
normalization factor which plays no role unless one wishes to compare different classes of 
models, so for the moment it can be dropped from the previous formula. One can say 
that in eq. (2.2) the first factor (the likelihood) is objective, while the second (the prior) 
contains our prejudices about how the probability is distributed a priori in the parameter 
space, given all our previous knowledge about the model. 

Ignoring the prior factor is not necessarily the most reasonable or "free of prejudices" 
attitude. Such procedure amounts to an implicit choice for the prior, namely a completely 
flat prior in the parameters. However, one needs some theoretical basis to establish, at 
least, the parameters whose prior can be reasonably taken as flat. 

Besides, note that a choice for the allowed ranges of the various parameters is necessary 
in order to make statistical statements. Often one is interested in showing the probability 
density of one (or several) of the initial parameters, say 9i, i = 1, N\, but not in the 
others, 9i, i = Ni + 1, ...,N. Then, one has to marginalize the latter, i.e. integrate in the 
parameter space: 

p(6i, i = l,...,iVi|data) = J d9 Nl+1 , ...,d6 N p{0 h i = 1, iV|data) . (2.3) 

This procedure is very useful and common to make predictions about the values of partic- 
ularly interesting parameters. Now, in order to perform the marginalization, we need an 
input for the prior functions and for the range of allowed values of the parameters, which 
determines the range of the definite integration (2.3). A choice for these ingredients is 
therefore inescapable in trying to make Bayesian LHC forecasts. 
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In consequence, an important indicator of the reliability of the results of a Bayesian 
analysis is the dependence of the results on the choice of the prior and on the ranges 
established for the unknown parameters, something that is addressed in the present paper 
(and very seldom done). 

2.2 Model-independent incorporation of naturalness criteria in the MSSM 

The main motivation of low-energy SUSY is the nice implementation of the electroweak 
(EW) breaking, since the EW scale (or, equivalently, the Z mass) does not suffer from 
enormous (quadratic) radiative corrections that occur in the SM. Actually, in the MSSM 
the EW breaking occurs naturally in a substantial part of the parameter space. This success 
is greatly due to the SUSY radiative contributions to the Higgs potential. Of course, in 
our analysis, the points of the parameter space that do not have a correct EW breaking 
are to be discarded, as usual. 

It is common lore that the parameters of the MSSM, {mo, Tni/%, A, B, ^}, should not be 
far from the experimental EW scale in order to avoid unnatural fine-tunings to obtain the 
correct size of the EW breaking. Previous Bayesian studies of the MSSM have attempted 
to incorporate this criterion by implementing some penalization of the fine-tuned regions 
or (more usually) by restricting the range of the soft parameters to < few TeV. This is 
more or less reasonable, but makes the results dependent on the actual ranges considered 
(although normally such test is not presented). 

On the other hand, since the naturalness arguments are deep down statistical argu- 
ments, one might expect that an effective penalization of fine-tunings should arise from 
the Bayesian analysis itself, with no need of introducing "naturalness priors" or restricting 
the soft terms to the low-energy scale. It was shown in ref. [12] that this is indeed the 
case (see also [13] for a previous observation in this sense). The key point is the following. 
Normally, the value of \i (up to a sign) is deduced from the value of the other parameters by 
the condition that the Higgs vacuum expectation value is the correct one, or equivalently 
the Z— mass is the experimental one. However, this procedure does not incorporate the 
fact that such value of \i may involve an extreme fine-tuning, in which case should become 
statistically penalized. This shortcoming can be surmounted by simply considering M^ xp 
as experimental data on a similar foot to the others, entering the total likelihood, C For 
the sake of simplicity let us approximate the likelihood associated to the Z mass as a Dirac 
delta, so 



where £ rcs t is the likelihood associated to all the physical observables, except Mz- Now, 
we can take advantage of this Dirac delta to marginalize the posterior in one of the initial 
parameters, e.g. /j,, performing a change of variable /x — > Mz- 



p(data|s, mo,mi/2,A, B,fj,) ~ S(Mz — M^ xp ) C 



rest i 



(2.4) 




(2.5) 
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where [iz is the value of [i that reproduces the experimental value of Mz for the given values 
of {s, mo, mx/2,A, B}, and p(s, tuq, mx/2,A, B, //) is the prior in the initial parameters (still 
undefined). Note that the Jacobian factor, can be written as where = 

is the conventional Barbieri-Giudice measure [14, 15] of the degree of fine-tuning. 



dlnMf 



<91n/x 

Thus, the above Jacobian factor has built-in the desired fine-tuning penalization, with no 
ad hoc assumptions. 

Actually this is enough to make the high-scale region of the parameter space, say 
soft terms > O(10) TeV, statistically insignificant; which allows in turn to consider a 
wide range for the soft parameters (up to the very Mx)- In consequence, the results of 
our analysis are essentially independent on the upper limits of the MSSM parameters, in 
contrast with previous studies. 

Beside the above [i — > Mz change of variable, it is highly advantageous for phenomeno- 
logical studies to trade the Yukawa couplings by the physical fermion masses (which are 
easily marginalized out) and the initial B— parameter by the derived tan/3 = (H u ) / (H c [) 
parameter. Consequently we should compute the whole Jacobian, J, of the transformation 

{n,y u B} -> {M Z) mt, tan/3} . (2.6) 

The initial prior (still undefined) is therefore always affected by this model-independent 
Jacobian factor, thus resulting in an effective prior in the new variables 

Pes(9i,mt,mo,m 1 / 2 ,A,tan(3) = J| M =^ p(g i ,yt,mo,m 1 /2, A, B, (jl = Hz) , (2.7) 

where we have already marginalized Mz using the associated likelihood ~ S(Mz — M^ xp ) 
(recall that [iz is the value of \i that reproduces the experimental Mz-) In eqs.(2.6, 2.7) 
we have made explicit the dependence just on the top Yukawa coupling and mass, but for 
other fermions goes the same. 

So, to prepare the scan in the new variables we need an explicit evaluation of the 
Jacobian factor. For the numerical analysis we have evaluated J using the Sof tSusy code 
[16] which implements the full one-loop contributions and leading two-loop terms to the 
tadpoles for the electroweak breaking conditions with parameters running at two-loops. 
This essentially corresponds to the next-to-leading log approximation; see ref. [17, 18, 19, 
20] for details. 

However, it is possible to give an analytical and quite accurate expression of J, and 
thus of the effective prior, by working with the tree-level potential with parameters running 
at one-loop (i.e. essentially the leading log approximation). At this level, the relation of 
{n, yt, B} with {Mz, rrtt , tan /3} is given by the minimization conditions and the expression 
for the (running) top mass, 

m^-tan 2 /3m^ M| 
^ " tan* 0-1 2" (2 - 8) 



Slow = 7^(m 2 Hl + m 2 H2 + 2^ ow ) (2-9) 
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yi ow 



V S/3 



(2.10) 



Here y denotes the top Yukawa coupling and the "low" subscript indicates that the quantity 
is evaluated at low scale (more precisely, at a representative supersymmetric mass, such as 
the geometric average of the stop masses). The soft masses m\ , m\ are also understood 
at low scale. For notational simplicity, we have dropped the subscript t from the Yukawa 
coupling. (Note that all these low-energy quantities contain an implicit dependence on the 
top Yukawa coupling through the corresponding RG equations.) Then 



J oc 



E 



tan 2 fj-1 



y\ ow tan f}(\ + tan 2 /3) ' 



(2.11) 



where Rn and E are defined through the renormalization group equations of \i and y, 



/'lo 



yE(Q\ ow ) 
1 + 6yF(Q 



low J 



(2.12) 



Here Q is the renormalization scale, F = jq^ h ElnQ, and E(Q) is a definite function that 
depends just on the gauge couplings [21] 1 . 

The important point is that the Jacobian (2.11) is a model-independent factor, valid 
for any MSSM, which must be multiplied by whatever prior is chosen for the initial 
parameters. It cannot be ignored. Hence, the effective prior in the standard variables 
{mt, mo, m.i/2, A, tan/3} reads 



p e f £ (m t ,m ,m 1 /2,A,tanf3) oc 



E 



tan 2 p-1 B 



low 



yi ow t(l + tan 2 /3) \l z 



p{m ,m l/2 ,A,B,^z), (2.13) 



where all the subjectivity lies in the p(mo,mi/2, A, B, fj,) piece, i.e. the prior in the initial 
parameters, for which we have still to make a choice (we will do it in the next subsection). 
This Jacobian (model-independent) prefactor contains the above-discussed penalization of 
fine-tuned regions. Thus, using it automatically incorporates the fine-tuning penalizations 
with no ad hoc assumptions. Note that, besides the penalization for large fj,, the Jacobian 
factor contains a penalization of large tan /3, reflecting the smaller statistical weight of this 
possibility. Actually, the implicit fine-tuning associated to a large tan (3 was already noted 
in ref. [22, 23], where it was estimated to be of order 1/ tan/3, in agreement with eq.(2.13). 



2.3 Priors in the initial parameters 

The choice of the prior in the initial parameters, {mQ,mi/2-, A, B, /j,} must reflect our knowl- 
edge about them, before consideration of the experimental data (to be included in the 
likelihood piece). Four of them, {m, M, A, B}, are soft SUSY-breaking parameters. They 
typically go like ~ F/A, where F is the SUSY breaking scale, which corresponds to the dom- 
inant VEV among the auxiliary fields in the SUSY breaking sector (it can be an F— term 



1 Besides an (irrelevant) numerical factor, the expression (2.11) contains a sin/3 factor which is cancelled 
assuming, as usual, a logarithmically flat prior for the Yukawa coupling, i.e. p(y) oc y~ ■ 
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or a D— term) and A is the messenger scale, associated to the interactions that transmit the 
breaking to the observable sector. Since the soft-breaking terms share a common origin it is 
logical to assume that their sizes are also similar. Of course, there are several contributions 
to a particular soft term, which depend on the details of the superpotential, the Kahler 
potential and the gauge kinetic function of the complete theory (see e.g. ref. [24]). So, it 
is reasonable to assume that a particular soft term can get any value (with essentially flat 
probability) of the order of the typical size of the soft terms or below it. The /i— parameter 
is not a soft term, but a parameter of the superpotential. However, it is desirable that its 
size is related (e.g. through the Giudice-Masiero mechanism [25]) to the SUSY breaking 
scale. Otherwise, one has to face the so-called /x— problem, i.e. why should be the size of 
ji similar to the soft terms', as is required for a correct electroweak breaking [see eq.(2.8)]. 
Thus, concerning the prior, we can consider fj, on a similar foot to the other soft terms. 

Let us denote Ms the typical size of the soft terms in the observable sector, Ms ~ F/A. 
Then, we define the ranges of variation of the initial parameters as 

- qM s < B < qM s 
-qM s < A < qM s 

< m < qM s 

< mi/2 < qM s 

< fj, < qM s (2.14) 

where q is an 0(1) factor. We have considered here the branch of positive \x. For the 
negative one we simply replace [i — > —fi. We can take q = 1 with no loss of generality, 
provided Ms is allowed to vary in the range < Ms < oo. In practice, to avoid divergences 
in the priors, we have to take a finite range for Ms, say 

M s < M s < M x , M^ ~ 10 GeV (2.15) 

Nevertheless, the values of the upper and lower limits of the Ms range are irrelevant. A 
small value for Ms (as the above lower limit) is in practice excluded by the LHC exclusion 
limits. A high value for Ms (as the above upper limit) is in practice statistically irrelevant 
thanks to the Jacobian prefactor in eq.(2.7) or eq.(2.13), i.e. by the naturalness of the EW 
breaking. So the results of the present approach are certainly independent of the ranges 
established for the parameters, something very satisfactory. 

Concerning the shape of the priors, as already stated, we find reasonable to assume 
(conveniently normalized) flat priors for the soft parameters inside the ranges (2.14), i.e. 

p{m ) = p(m 1/2 ) = p(n) = — , p(A) = p(B) = (2.16) 

The next step is to choose a prior for Ms- Admittedly, this is the less objective part 
of the statistical analysis, but, as discussed above, one cannot simply ignore the prior; 
it is an unavoidable piece when one wants to determine the probability distribution in 
the parameter space. The most conservative attitude is to use two different, though still 



- 8 - 



reasonable, priors, and then compare the results. This gives a fair measure of the prior- 
dependence. 

For that matter we have considered two somehow standard types of prior: flat and 
logarithmic. In a flat (logarithmic) prior one assumes that, in principle, the typical size 
(order of magnitude) of the soft terms can be anything, say from Mg up to Mx, with 
equal probability. In our opinion, a logarithmic prior is probably the most reasonable 
option, since it amounts to consider all the possible magnitudes of the SUSY breaking in 
the observable sector on the same foot (this occurs e.g. in conventional SUSY breaking by 
gaugino condensation in a hidden sector). However, we will consider log and flat priors at 
the same level throughout the paper, in order to compare the results and thus evaluate the 
prior-dependence. 

Finally, we can marginalize Mg, which thus disappears completely from the subsequent 
analysis, leaving a prior which depends just on the {mo, mi; 2) A, B, fj,} parameters 2 . Using 
eqs.(2.16), one gets, for the logarithmic case, 

■Mx 



1 



p(m ,m 1/2 ,A,B,fi) oc / ^ dM s 



max{m ,m 1/2 ,\A\,\B\,ij,,M°} ^ S 
1 



(2.17) 



[max{m ,m 1/2 , \A\, \B\, [i, M®}] 5 

which is the prior to be plugged in eq.(2.7) (or in the approximated expression (2.13)) to 
get the effective prior in the scan parameters. For the flat case, the analogous computation 
gives 

P (,„„, m 1/2 , A, B, „) oc [max{m0imi/2i ^ UBU , M „ )]4 . (2-18) 

The fact that eqs.(2.17, 2.18) correspond, respectively, to logarithmic and flat priors in 
the size of the soft terms can be verified by marginalizing all the soft parameters but 
one, say m^, before incorporating any experimental information. Then one does obtain 
logarithmic and flat behaviors, respectively, as is clear from the exponents in the above 
denominators (for a more throughout discussion see ref. [27]). 



3. Physical observables and computational methods used in the analysis 

A first group of physical observables used in the analysis are the SM-like parameters, e.g. 
masses and couplings of SM particles. Some of them, in particular the top and bottom 
masses and the electromagnetic and strong coupling constants, have a strong influence in 
the analysis. We have treated them as nuisance variables (marginalizing them) using the 
central values and uncertainties given in Table 1 

Furthermore, the full set of experimental data used to compute the likelihood function 
is listed in Table 2. The full likelihood function is the product of the individual likelihoods 

2 This procedure is a "hierarchical Bayesian technique", first used in ref. [26], but using complicated 
functions that were not possible to integrate analytically. 
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SM nuisance parameters 




Gaussian prior Range scanned 


ref. 


M t [GeV] 


173.2 ±0.9 (167.0,178.2) 


[28] 


m b {m b ) Ms [GeV] 


4.20 ±0.07 (3.92,4.48) 


[29] 


[OemiMz)^ 3 ]- 1 


127.955 ±0.030 (127.835,128.075) 


[29] 


a s (M z ) Ms 


0.1176 ± 0.0020 (0.1096, 0.1256) 


[30] 



Table 1: Nuisance parameters adopted in the scan of the CMSSM parameter space, indicating the mean 
and standard deviation adopted for the Gaussian prior on each of them, as well as the range covered in the 
scan. 

associated to each piece of experimental data. As discussed in the next section, not all 
those observables are included in the computation of the full likelihood all the time. In 
particular, we plug and unplug the dark matter constraints to check its impact. For the 
quantities for which positive measurements have been made (as listed in the upper part of 
Table 2), we assume a Gaussian likelihood function with a variance given by the sum of 
the theoretical and experimental variances, as motivated by eq. (3.3) in ref. [31]. For the 
observables for which only lower or upper limits are available (as listed in the bottom part 
of Table 2) we use a smoothed-out version of the likelihood function that accounts for the 
theoretical error in the computation of the observable, see eq. (3.5) and fig. 1 in ref. [31]. 

Let us comment briefly on some of the observables and the theoretical method used to 
compute the theoretical prediction for them and the corresponding likelihood. 

To calculate the MSSM spectrum we have used SoftSusy [16], where SUSY masses are 
computed at full one-loop level and the Higgs sector includes twodoop leading corrections 
[44]. We discard points suffering from unphysicalities: no self- consistent solutions to the 
RGEs, no EW breaking and tachyonic states. In our treatment of the radiative corrections 
to the electroweak observables My/ and sin 2 9 e s we include full twodoop and known higher 
order SM corrections as computed in ref. [45, 46], as well as gluonic twodoop MSSM 
corrections obtained in [47]. 

The most important piece of experimental information is by far the Higgs mass, rrih, 
for which we use the value obtained from the combination of the ATLAS and CMS experi- 
mental results, expressed in eq.(1.2). At each point of the supersymmetric parameter-space 
we evaluate the theoretical prediction for rrih through SoftSusy, which has the two-loop 
computation incorporated. The estimated theoretical uncertainty is ±2 GeV. 

Concerning B-physics, the branching ratio for the B — > X s j decay (the most important 
one) has been computed with the numerical code SusyBSG [48] using the full NLO QCD 
contributions, including the two-loop calculation of the gluino contributions presented in 
[49] and the results of [50] for the remaining non-QCD tan /3-enhanced contributions. The 
supersymmetric contributions to b — > 57 grow with decreasing masses of the supersymmet- 
ric particles and with increasing tan f3. The other B(D)-physics observables summarized in 
Table 2 have been computed with the code Superlso (for details on the computation of 
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Observable 


Mean value 


Uncertainties 


ref. 






a (exper.) 


r (theor.) 




M w [GeV] 


80.399 


0.023 


0.015 


[32] 


sin 2 6 P f f 


0.23153 


0001 6 
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m ,m 1/2 


LHC exclusion limits, 


see text 


[41] 


rriA, tan/3 


LHC exclusion limits, 


see text 


[42] 


SI 


XENON100 2012 limits (224.6 


x 34 Kg days) 


[43] 



Table 2: Summary of the observables used for the computation of the likelihoo d funct ion For each 
quantity we use a likelihood function with mean /i and standard deviation s — V 'a 2 + r 2 , where a is 
the experimental uncertainty and r represents our estimate of the theoretical uncertainty. Lower part: 
Observables for which only limits currently exist. The explicit form of the likelihood function is given in 
ref. [31], including in particular a smearing out of experimental errors and limits to include an appropriate 
theoretical uncertainty in the observables. 

the observables see [51] and references therein). Among them is the BR{B S —> u + u~) for 
which we include the new LHCb measurement (3.2^2) x 10~ 9 , derived from a combined 
analysis of 1 fb _1 data at y/s = 7 TeV collision energy and 1.1 ftr 1 data at *Js = 8 TeV 
collision energy [38]. We implement this constraint as a Gaussian distribution with a con- 
servative experimental error of a = 1.5 x 10 -9 , and a 10% theoretical error. Both codes 
have been integrated into the forthcoming SuperBayes-v2 . code [52]. 

We also include recent data from the XENON100 direct detection experiment, obtained 
with an effective volume of 34 kg and an exposure of 224.6 days between February 2011 and 
March 2012 [43]. Two candidate scattering events were detected within the signal region, 
with an expected number background of&=1.0±0.2 events, resulting in exclusion limits 
in the (m^o, a^_ p ) plane. The corresponding likelihood function is described in detail in 
ref. [11], which we refer the reader to for further details. 

In addition, we incorporate the most recent constraints from SUSY searches at the 
LHC. We follow the analysis presented by the ATLAS collaboration in ref. [53], since it 
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is the one associated to the more stringent constraints: a data sample corresponding to a 
luminosity of 5.8 fb _1 for a center-of-mass energy of \/s = 8 TeV in the 0-lepton channel. 
As no significant deviations have been so far found with respect to SM predictions, this 
has been used to derive the constraints on the (mo, m^) plane [41]. We note that while 
this 95% exclusion line has been obtained assuming values of tan (3 = 10 and Aq = GeV, 
the final result is not expected to change much for different values of these two parameters 
since the production and decay of squarks and gluinos are quite insensitive to the values of 
tan /3 and Aq [54] . Besides, we have applied this limit to the NUHM since the SUSY signal 
for gluino and squarks production does not depend significantly on other parameters than 
mo and m^. Additionally we have included exclusion limit in the (m^jtan/?) plane from 
a CMS search for the decay of neutral Higgs bosons into final states containing two muons 
and missing Et, based on 4.5 fb _1 integrated luminosity of data collected at yfs = 7 TeV 
collision energy [42]. We include those exclusion limits by assigning a likelihood equal to 
zero to samples falling below the exclusion line. 

As a scanning algorithm we adopt the MultiNest [55] algorithm as implemented in 
SuperBayeS-v2 . 0. It is based on the framework of Nested Sampling, invented by Skilling 
[56, 57]. MultiNest has been developed in such a way as to be an extremely efficient sampler 
even for likelihood functions defined over a parameter space of large dimensionality with a 
very complex structure as it is the case of the CMSSM. The main purpose of the Multinest 
is the computation of the Bayesian evidence and its uncertainty but it produces posterior 
inferences as a by-product. For the marginalization procedure we have used the above- 
discussed ranges for our priors. Besides, we have considered 2 < tan/3 < 62, although the 
precise limits of this range are irrelevant. A too low value of tan f3 is excluded by the Higgs 
mass value, unless the supersymmetric masses are extremely high, which is drastically 
disfavoured by the naturalness of the EW breaking through the above-discussed Jacobian 
factor. (More precisely, tan/3 = 2 requires stop masses above 10 3 TeV [58, 59]). On the 
other hand, the region of very large tan f3 is also strongly suppressed by the EW mechanism, 
as is clear by the prefactor in eq.(2.13). 

4. Results for the CMSSM 

We have organized the presentation of the results in following way. 

In subsect. 4.1 we show the plots of the probability distribution function (in one and 
two dimensions) using only the most robust group of experimental data, namely LHC 
exclusion limits, Higgs mass, and EW- and B(D)-physics observables. At this stage we 
leave aside dark matter constraints since supersymmetric dark matter could be unestable 
with a tiny amount of R-parity violation (without modifying the rest of the analysis). In 
that case the dark matter constraints would be irrelevant. In subsect. 4.2 we include the 
dark matter constraints in two different fashions: either requiring that the supersymmetric 
dark matter is consistent with actual dark matter observed abundance or less than it. Note 
that the latter, more conservative, case could be realistic if there is an additional source of 
dark matter (such as axions). In subsect. 4.3 we discuss the g-2 issue. Along this section 
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we present the results just for the CMSSM scenario. The results for the NUHM one are 
presented in section 6 following a similar scheme. 



4.1 Results with no dark matter constraints 

At this stage we have used all the observables, except those related to dark matter. 

In fig. 1 we present the probability distribution function (posterior) for the four CMSSM 
parameters, mo, nti/2) A tan/3, and the physical gluino and (first and second generation) 
squark masses. For each case, the non-plotted parameters have been marginalized. The 
red(cyan) lines show the results for flat (logarithmic) priors All the plots in fig. 1 correspond 
to positive //. The results for negative ^ are almost identical, so we do not show them to 
avoid proliferation of figures. 

Fig. 2 shows the probability distribution functions for flat and log priors when all 
parameters except two are marginalized. This gives rise to the two-dimensional plots in 
the mo — m\j2 and mq — m g planes. The present LHC limits (at 5.8 fb _1 )[41] and the 
expected ones in the future (3000 fb~ 1 )[60] are shown for reference. 

Clearly, the bulk of the probability lies in the high-energy region which is not directly 
accessible to the LHC at present. The rather high Higgs mass reported by ATLAS and CMS 
is the main responsible of this effect. As it is well-known, in the MSSM the tree-level Higgs- 
mass in the MSSM is bounded by the mass of the Z-boson, so large radiative corrections 
are needed in order to reconcile theory and experiment. An approximate analytic formula 
for ffiji [61, 62, 63, 64, 65, 66, 20, 67], useful for discussions, reads 



to^ Mf cos 2 2/3 + 



log M susy , *t ( 1 X± 

m t JW SUSY V 1ZJW SUSY 



+ • • • (4.1) 



where m t is the top running mass, Msusy represents a certain average of the stop masses 
and Xt = At — fi cot /3 [65, 20, 67]. The first term of (4.1) is the tree-level Higgs-mass and 
the second are the dominant radiative and threshold contributions. Note that the radiative 
corrections grow logarithmically with the stop masses while the threshold correction has a 
maximum for Xt = ±\/6MgusY- 

The important point is that, in order to achieve to^ = 125 — 126 GeV, one typically 
needs (besides non-small tan/3) stop masses > 3 TeV, unless Xt happens to be close to the 
above-mentioned maximizing value (we will comment on the latter possibility below) . This 
is what the probability distributions of Figs. 1, 2 show. It should be kept in mind that in 
the CMSSM all squark masses come from a universal value at Mx, so it is not possible to 
have large stop masses (required to reproduce the value of m/,) and lighter masses for the 
other squarks. In a general MSSM Figs. 1, 2 would approximately represent the constraints 
on the third-generation squark masses and trilinear coupling, and on the gluino mass (at 
Mx and at low-energy). 

Of course, these kinds of results are known and expected; but the present Bayesian 
analysis allows to quantify the effect, showing the new probability distribution functions 
of the various parameters. 
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Figure 1: ID marginalized posterior probability distribution of the CMSSM parameters: 
mo, Tii/2, A, tan/3 and the physical squark (of the first two generations) and gluino masses. The 
rcd(cyan) corresponds to flat (logarithmic) priors. 



The results of Figs. 1, 2 show clearly that, with the present information, the probability 
of detecting the CMSSM at present at the LHC (even if it is the true model) are < 5%, 
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Figure 2: 2D marginalized posterior probability distribution of the CMSSM in the m 1 / 2 — mo and 
the rrig — rriq planes. Left panels show the flat priors case and right panels the logarithmic one. 
Both for /i > 0. The inner and outer contours enclose respective 68% and 95% joint regions. 



at least with the current techniques to search for physics beyond the SM in LHC data. 
In the long-term future the probability is substantially higher, but still the chances of not 
discovering it are larger. To put this in optimistic/wishful-thinking terms, since rrih ~ 126 
GeV, if one believes in the CMSSM one should not be worried about the no-detection of 
SUSY at LHC up to now! But its detection in the future will be challenging too. This is 
the main conclusion of the Bayesian analysis at this stage. 

On the other hand, it should be kept in mind that in a Bayesian analysis one assumes 
that the scenario studied (in this case the CMSSM) is the true one, and then one determines 
the probability distribution in the parameter space. It does not give a direct measure of 
how likely is the scenario itself, unless one compares it with an alternative scenario using 
evidence-comparison techniques. We will come back to this point in the discussion of sect. 
6 below. But it is already clear that the experimental results on the Higgs mass have 
pushed the probability to a region which was statistically disfavored (fine-tuned) before 
the discovery, which sheds shadows on the naturalness of the CMSSM. In sect. 6 we will 
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make this point more quantitative. 

Let us comment now on other interesting aspects of the results shown in Figs. 1, 2. As 
already mentioned, the results for negative [i are almost the same. The experimental Higgs 
mass requires heavy supersymmetric particles, which makes SUSY decoupled for many 
purposes and makes all the sign fi effects in observables (e.g. electroweak and B— physics 
observables) rather irrelevant. An exception to this would be the g — 2 data, which have 
not been included in the analysis for the reasons explained in subject. 4.3. On the other 
hand it is quite remarkable the stability of the results against changes in the priors. Al- 
though the log priors prefer logically slightly lower values of the supersymmetric masses, 
the difference is not significant. This is quite satisfactory, showing that the experimental 
data are already able to determine pretty unambiguously the preferred supersymmetric re- 
gions in the parameter space. Still, how ever, the acceptable ranges for the supersymmetric 
parameters are quite wide, but systematically beyond 1 TeV. 

The probability distribution for tan /3, with a maximum around tan /3 ~ 6, is originated 
by two opposed effects. On the one hand, the large experimental Higgs mass favors large 
tan/3, since the tree-level contribution to nih is proportional to cos 2 2/3 and therefore grows 
with tan/3; although the effect is negligible beyond tan/3 ~ 8. On the other hand the 
Jacobian factor (2.11) penalizes large tan/3, due to the implicit fine-tuning among the 
initial parameters required to obtain it. 

The Aq parameter shows two wide peaks for positive and negative values, see fig. 1. 
This effect is in fact due to the logarithmic scale used for the representation. In a linear 
scale, the probability presents just one peak in the central region. On the other hand, it 
can be noticed that the preferred positive values of Ao tend to be larger than the negative 
ones. This comes from the fact that the RG equations drive A towards negative values, 
which increases (decreases) its absolute value if its initial sign is negative (positive). Thus 
a too large initial negative value of Aq gives too large values of A at low energy, leading to 
the appearance of tachyons. 

The discussion of the A— parameter raises another important point. As mentioned 
above, the value of rrih is maximized for Xt = ±\/6Msusy- In that case, lower masses 
for the stops are acceptable, which in principle would make SUSY less fine-tuned. This 
possibility is certainly interesting. However, it represents itself another kind of fine-tuning 
unless one has some theoretical reason to expect Xt in the favorable range. Incidentally, 
note that for moderately large tan /3 the value of Xt is essentially given by At, so one would 
essentially need a value of At in the favorable range. In any case, this would amount to go 
beyond the MSSM, e.g. invoking a particular SUSY breaking and mediation mechanism 
that could lead to this result for At 3 . If we do not rely on a mechanism of that kind, the 
question is whether it pays or not (in statistical weight) to have Xt around the maximal 
value. This question is answered by the Bayesian analysis. If the answer is positive, 
the value of Xt should be concentrated around that privileged value. Fig. 3 shows the 
probability distribution function in the plane — (Xt/Ms) (all the remaining parameters 

3 See ref.[68] for a paper arguing that this can actually occur in a class of superstring scenarios. 
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Figure 3: 2D marginalized posterior probability distribution of the CMSSM in the mh-X t /Ms 
plane, using log priors for /i > 0. The inner and outer contours enclose respective 68% and 95% 
joint regions. 



are marginalized). Clearly, the probably does not show any preference for that value. In 
other words, requiring that value for Xf implies a more severe fine-tuning than living with 
rather large squark masses. 

4.2 Incorporation of dark matter 

There are different astrophysical and cosmological observations that offer impressive ev- 
idence of the existence of dark matter (DM) in the Universe (see Table 2 for a recent 
determination of ^dm)- On the other hand, the consistency with the observed large struc- 
ture of the Universe favours cold dark matter (CDM), i.e. non-relativistic matter at the 
beginning of galaxy formation. This leads to the hypothesis of a weakly interacting massive 
particle (WIMP) as the component of CDM. 

Supersymmetry offers a good candidate for such a WIMP, namely the LSP, which is 
stable in the standard (R-parity conserving) SUSY formulations (for a review see [69]). 
Although, depending on the models, there are several possibilities for the SUSY WIMP, 
the most popular and natural candidate is the lightest neutralino (from here on we will refer 
it as the neutralino), x?> which is the LSP in most of the CMSSM parameter space. One 
could assume that the LSP makes up the whole content of DM required by the "standard" 
cosmology. However, it could happen that the early thermal production of WIMPs is 
diluted by electroweak inflation. Admittedly, the latter is not a most natural or popular 
scenario of inflation, but mechanisms for it have been explored [70]. Alternatively, the LSP 
could be unstable assuming tiny violations of R-parity, see e.g. [71]. In these cases the 
observed dark matter should be provided by other candidate, e.g. an axion. 

Motivated for this, we discuss next, separately, the cases where the DM is made only 
of LSPs and where it is partially made of LSPs. 
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Single- Component DM 



In this scenario we assume a Gaussian likelihood on Q x h 2 (see Table 2). Besides, we 
apply the constraints from DM-induced recoil events in the XENON 100 experiment, as 
described above. For the present analysis we neglect the effect of the astrophysical and 
hadronic uncertainties entering into the determination of the spectrum of recoil events. 
Namely, we fix the astrophysical parameters that describe the density and velocity distri- 
bution of DM particles at the commonly adopted benchmark values: local CDM density 
A3,cdm = 0.4GeVcm -3 , circular velocity vo = 235 km s _1 and escape velocity v esc = 550 
km s _1 (see, e.g., [72] and references therein for a recent discussion of the astrophysical 
uncertainties on these quantities) . For the contribution of the light quarks to the nucleon 
form factors we have adopted the values fx u = 0.02698, frd = 0.03906 and fxs = 0.36 
[73]. Incorporating the associated uncertainties of both the astrophysical and nuclear form 
factors parameters to the analysis as nuisance parameters (i.e. marginalizing) would only 
mildly change the numerical results presented here; the main conclusions would remain 
unchanged, in agreement with the findings of [74]. 

In fig. 4 we plot the ID pdfs of various supersymmetric parameters and the physical 
gluino and squark masses. The cyan-continuous (blue-dashed) line corresponds to [i > 
(/j, < 0). This is complemented with the 2D pdfs of fig. 5. In all cases the non-displayed 
parameters have been marginalized. The results correspond to log priors (for flat priors 
they are very similar). The first thing to notice is that, due to the DM constraints, the 
pdfs of m\ji and rriQ (an thus of nig, nig) have moved further to large values. Indeed, more 
than the 95% of the posterior probability is out of the LHC reach, even with an energy 
of 14 TeV and a integrated luminosity of 3000 fb _1 as projected for the HL-LHC (High 
Luminosity) phase. [60]. This conclusion holds for the flat priors case. The reason is the 
following. There are in fact two additional regions where the required DM abundance can 
be reproduced: the stau-coannihilation and the A— funnel regions. However, the first one is 
disfavored, as it corresponds to rather small soft terms, difficult to reconcile with the Higgs 
mass. The second one requires intermediate/large tan /3 values, which are penalized by the 
Jacobian prefactor in the r.h.s. of eq.(2.13). In consequence, one is essentially left with the 
Focus-Point region (actually the other two regions survive but with very tiny statistical 
weight, as discussed in sect. 7). 

In Appendix A we describe some characteristics of the Focus-Point in the CMSSM and 
the NUHM. Let us now mention some of them, which will be useful to understand the shape 
of Figs. 4, 5. The Focus- Point region corresponds to heavy scalars (mo), less heavy gauginos 
( m i/2) an d tan/3 > 8. In that region the value of rrifj at low energy is quite insensitive to 
its value at high-energy, this is the reason why it can reach large values with no additional 
fine-tuning. For too-large mo, however, the electroweak breaking solution gets lost. This 
occurs when the required value of /j, 2 gets negative. For slightly smaller values of mo the 
value of /u is small, and this is welcome for fine-tuning reasons, since the Jacobian prefactor 
in eq.(2.13) gets smaller. Hence, this is a theoretically interesting region by itself. But, 
besides, it is also appropriate to correctly reproduce the DM abundance. The reason is that 
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Figure 4: As in fig. 1 but with an additional constraint from DM (WMAP 7-years and XENON100) 
considering that the lightest supersymmetric particle (LSP) is the soley DM component. Here the 
cyan line corresponds to a positive fi whereas the blue dashed to a negative /x. 



if ji is smaller than (or similar to) Txi\/2i the light est neutralino is essentially a Higgsino (or 
a mixture of Higgsino and bino). That sizable fraction of Higgsino makes the neutralinos to 
coannihilate efficiently with the lightest chargino and the second lightest neutralino, which 
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Figure 5: As in fig. 2 but with an additional constraint from DM (WMAP 7-years and XENON100) 
considering that the LSP is the soley DM component. 



may lead to the correct DM abundance. Actually, if the lightest neutralino is essentially a 
pure Higgsino and its mass is too small, the annihilation is too efficient. Then one needs 
a sizable bino component to reduce the annihilation so that the correct DM abundance 
is reproduced. Therefore the largest acceptable neutralino mass corresponds to the pure 
Higgsino case, and is M x ~ 1 TeV [75] (in that limiting case the neutralino, chargino, and 
second neutralino masses are all ~ Below that limit, the neutralino mass can vary in a 
wide range, provided it contains the suitable bino fraction. Indeed, this was the situation 
prior to the new XENON100 data. 

The point is that XENON is sensitive to the spin-independent cross-section of O(100) 
GeV neutralinos with the nucleons, if the former contain both a Higgsino and bino com- 
ponents. This occurs through diagrams exchanging a Higgs state. In this way XENON100 
has been able to exclude all the above-mentioned Focus-Point region, whenever there is 
a non-negligible Higgsino-bino mixture. In essence, only the pure Higgsino possibility re- 
mains alive [5]. In consequence the value of m x and \i are constrained around 1 TeV. 
In addition, the bino mass cannot be smaller than that, otherwise the lightest neutralino 
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Figure 6: 2D marginalized posterior probability distribution of the spin-independent scattering 
cross section of a neutralino off a proton versus neutralino mass being this the soley DM component 
for log priors, [i > (left panel) and fi < (right panel). 



would get a sizable bino component, which is excluded by XENON. This is the reason of 
the lower limit, rrii/2 > 2 TeV, clearly visible in Figs. 4, 5 (the bino mass is approximately 
half the initial gaugino mass at Mx)- 

Although the prospects of LHC to discover or discard this scenario are clearly pes- 
simistic, future XENON measurements could do the job. In fig. 6 we display the 2D joint 
posterior pdf of the spin-independent scattering cross section of a neutralino off a pro- 
ton versus the neutralino mass for log priors and both signs of fj,. As expected from the 
previous discussion, the posterior pdf of the neutralino masses is confined around 1 TeV. 
In contrast, the spin-independent scattering cross section covers several orders of magni- 
tude. The present XENON100 and future XENON1T limits are also indicated in the plots. 
Clearly, most of the significant region favoured at the 95 % level will be within reach of 
the next generation of direct detection experiments, which are certainly good news for the 
CMSSM. 

Let us finally notice the fact, that, contrary to what happened before including DM 
constraints, the pdf in Aq shows a clear trend towards positive values. This has to do with 
the Focus-point regime. If Aq is negative, its absolute value at low-energy increases. This 
increase is translates to the running of m 2 H ^ , which becomes more negative at low-energy. 
As a consequence, the required value of fj? to achieve successful electroweak breaking in- 
creases. This implies a larger fine-tuning and makes more difficult the implementation of 
the Focus-Point solution to the DM abundance. For positive Aq, the behavior is the op- 
posite. Without the DM constraints the Focus- Point area was already inside the preferred 
region, but it was not the only one, as can be seen from fig. 2, thus the two signs of Aq 
had approximately the same statistical weight. 

Multi- Component DM 
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Relaxing the requirement that all of the DM is made of LSPs implies that the WMAP 
measurement is just an upper bound on its abundance. At it is shown in the Appendix of 
ref. [76], in this case, the correct effective likelihood is given by the expression 

AvMAp(ttLSp/i 2 ) =Co e-^-^x^dx, (4.2) 

where Co is an irrelevant normalization constant, r* = /^wmap/cwmap and Q^sph 2 is the 
predicted relic density of the LSP as a function of the CMSSM and SM parameters. 

When the LSP is not the only constituent of DM, the rate of events in a direct-detection 
experiment is smaller, since it is proportional to the local density of the LSP, p x , which is 
now smaller than the total local DM density, pdm- The suppression is given by the factor 
£ = Px/pdm- Following ref. [77], we assume that ratio of local LSP and total DM densities 
is equal to that for the cosmic abundances, thus £ = p x /pdm = ^x/^dm- For CIbm we 
adopt the central value of the WMAP determination see Table 2, while for pdm we adopt, 
following ref. [72], the value 0.4GeVcm -3 . This allows to evaluate £ for each point in the 
CMSSM parameter space. 

The effect of allowing a multi-component DM in the ID posterior pdfs is displayed 
in Figs. 7, 8, which should be compared with the previous Figs. 4, 5, when the LSP was 
assumed to be the only component of DM. The most notable effect is that now large areas 
with lighter are recovered. This is easy to understand. As mentioned above, if the 

neutralino is pure Higgsino, for m x < 1 TeV, the LSP abundance is less than the observed 
DM abundance. This does not represent now any problem, as one assumes that the DM 
has other components. Hence, smaller values of m x , or equivalently of (j,, are now rescued. 
The condition for the Higgsino to be the LSP is that the bino mass is not smaller than 
\x. Consequently, m^j (and thus m g ) can now be smaller, as reflected in Figs. 7, 8. On 
the other hand, it is now possible for m^ 2 to go even below that limit. The reason is 
that, although in that case the neutralino becomes a mixture of Higgsino and bino, which 
is potentially dangerous from XENON limits, the rescaling factor £ reduces the predicted 
signal of recoil events. The actual value of £ depends on the neutralino mass and the 
composition of the mixture. 

The effect of all this is to bring m^ 2 to lower values, leading to gluino and squark 
masses of 0(1) TeV within the LHC reach. This can be seen in the lower panel of fig. 7. 
Therefore we conclude that, regarding LHC searches, there are much better detection 
prospects in the case of multi-component DM than in the single-component DM scenario. 

Finally, in fig. 9 we show the probability in the (m^o,^ 1 ) plane. The shape of the 
distribution is easy to understand from the above discussion. Neutralino masses of 0(100) 
GeV become now allowed at the 95 % Bayesian credibility level, though the upper bound 
remains at of ~ 1 TeV. In addition, values of cr^ 1 > 10 -8 pb are now within the 95 % 
level due to the suppresion caused by the £ factor in predicted recoil events signal. The 
XENON1T expected 90 % exclusion limit contour is also displayed, showing that the 
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Figure 7: As in fig. 4 but considering that the LSP is not the soley DM component. 



region favored by the data at the 95 % level can be fully probed by future tone-scale direct 
detection experiments. 
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Figure 8: As in fig. 5 but considering that the LSP is not the soley DM component. 
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Figure 9: As in fig. 6 but considering that the LSP is not the soley DM component. 



4.3 g-2 

The magnetic anomaly of the muon, = ^(g — 2)^ has been a classical and powerful test 
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for new physics. The theoretical determination of depends crucially on the evaluation 
of the contribution coming from the hadronic vacuum-polarization diagram. This can be 
expressed in terms of the total hadronic cross section e + e~ — > had or, alternatively, in 
terms of the hadronic r— decay data, which are theoretically related to the e + e~ hadronic 
cross section. Both approaches lead to different results, although in the last times they 
have notably converged. On the other hand, in either case the results are in remarkable 
disagreement with the experimental determination. More precisely, defining 5a^ = 5 e ^ v a^ — 
6 dp, the recent results are Sa^ = (25.5 ± 8.0) x 10~ 10 (using e + e~ data) and = 
(15.7 ± 8.2) x 1CT 10 (using r— decay data) [78, 79]. This corresponds to a discrepancy of 
3.2cr and 1.9a respectively. 

This potential signal of new physics may be attributed to supersymmetric contributions 
to the muon anomaly, 6 a^, which arise mainly from 1-loop diagrams with chargino 
or neutralino exchange [80, 81, 82, 83, 84]. This contribution i ncreases with increasing 
tan/3 and decreasing supersymmetric masses. The problem is that the large Higgs mass 
revealed by the experiment, well above Mz, requires large supersymmetric masses, which 
is hardly conciliable with a sizeable 5 MSSM a M . This is the reason why we have not used 
the g — 2 result in the computation of the likelihood in our analysis. Besides, the g — 2 
computation has been rather controversial in the last years and probably the situation is 
not clear enough yet. 

The Bayesian approach allows for a quantification of the discrepancy between the 
observed Higgs mass and the assumption that the supersymmetric contributions are re- 
sponsible for Sa^. We briefly describe the procedure (for more details see ref.[85, 86, 87]), 
which will also be used in sect. 6 in other context. The idea is to separate the complete set 
of data in two subsets: 



Here & represents the subset of observations, whose compatibity with the rest of the 
observations, D, we want to test. In our case, @ is the experimental value of a^, whereas 
D denotes the remaining experimental data, including the experimental Higgs mass. Then 
one constructs the quantity p(S>\D), i.e. the probability of getting a certain result when 
measuring 3> = a^, given the known values of the remaining observables, 



Here p(£), D) is probability of measuring both sets of data at the same time, while p(D) 
is the analogue, but just for the D subset. These probabilities are given by the integral 
of the likelihood times the prior in the parameter space, and therefore coincide with the 
normalization factor appearing in denominator of eq.(2.2), sometimes called the evidence. 



{data} = D} 



(4.3) 



p{9\D) 




(4.4) 



p{D) 



E.g. 



p(@, D) = p(data) 




(4.5) 
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Now, the consistency of £^ obs (i.e. the actually measured muon anomaly) with the rest of 
data, D, can be tested by comparing the probability p(2t ohs \D) with the one obtained for 
the value of @ that maximizes such probability, say f^ max (assuming the same reported 
error at the new central value). ^ max is somehow the "prediction" of the model for the 
muon anomaly, given the rest of experimental data. The quantity 

X(2*>\D) = P{ ^ D) = Pi ^ D) (4 6) 

is a global likelihood ratio (integrated to all possible values of the parameters of the model) . 
Note that the p(D) factor cancels out in the expression of i?(f^ obs |D), which is simply 
given by the ratio of the global evidences. This was called the ££— test in ref.[87]. If this 
likelihhod has a Gaussian-like profile, the quantity — 2 lnif (ii? obs |L>) can be interpreted as 
g— deviation from the theoretical prediction. 

Obviously, for very large m^, the supersymmetric particles must be very heavy. In that 
limit SUSY decouples, so ^ max should approach the SM prediction, i.e. -2 In if (S> ohs \D) -> 
3.2 2 (asuming the theoretical calculation of a M based on e + e _ data); which is indeed the 
case [88]. Now, one can compute if for any value of m^, in particular for = 125.9 ±0.4 
GeV, as given by the experiment. Our calculation indicates that, essentially, for this value 
of the Higgs mass, if has already reached its maximum (decoupling limit) value. In other 
words, the tension between the experimental Higgs mass and the experimental a M is already 
as large as for the SM (and thus larger than 3a). This means that in the regions of the 
CMSSM where the prediction for is closer to the experiment than the SM prediction are 
not favoured since they present an onerous disagreement with t he observed Higgs mass. 
In some tiny regions this price can be somewhat alleviated, but they are statistically quite 
irrelevant, which reflects the cost of the underlying fine-tuning. 

5. NUHM 

Let us now (more briefly) describe the results when the MSSM model considered is the Non 
Universal Higgs Masses one (NUHM). Most of the discussion (and many of the results) are 
in fact analogous and it is not necessary to repeat them here. Hence, we will focus on the 
new features and their origin. 

A slight difference has to do with the choice of the priors. Actually, the priors we have 
used for the NUHM are the same as the ones described in eq. (2.14) - (2.18), but now 
we have to include the common mass for the Higgses at Mx, i-e. mn as an additional 
independent parameter. Then, the log prior reads, 

p(m ,m ff ,m 1/2 , A, B,fj, oc ^ max j m ^ m ^ m 1/2 , \A\, \B\, //, Mg}] 6 ' 
and the flat prior, 

p(m ,m H ,m 1/2 ,A,B,n) oc ^ToTu- ( 5 - 2 ) 

[max{m ,m H ,m 1 /2, \A\, \B\, fi, M£}J 4 
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Figs. 10, 11 show the results when DM constraints are not included in the analysis. They are 
analogous to Figs. 1, 2 for the CMSSM. Fig. 10 shows the posterior probability distribution 
function for mo, run, m\/2i an d tan/3, for log and flat priors and \i > 0. The results for 
fi < are quite similar and are not shown. As in the CMSSM, the preferred region is in 
the high-energy region, mainly due to the large value of the Higgs mass. From the mo and 
run pdfs in fig. 10 or from the two-dimensional pdf in the m# — mo plane in fig. 11 it is 
clear that the NUHM naturally prefers mo < win- The reason is that, as discussed in the 
Appendix, the separation of the mo and ttih values can be used to decrease the value of mo 
at which [i 2 becomes very small (and eventually negative). Recall that in the CMSSM this 
occurred for large mo, and corresponded to the region of high mo in fig. 2. Actually, after 
including DM constraints, that "Focus-Point" region was the only surviving one. Now, it 
is possible (and pays) to live with smaller mo, which statistically is advantageous. This is 
the origin of the low— mo region in Figs. 10, 11. 

Finally, Figs. 12, 13 show the results after including DM constraints (assuming that 
the DM is purely LSP). The results are quite similar to those for the CMSSM, but shifted 
to the low-mo region for the above-commented reasons. The same happens if we allow 
multi-component DM (we do not show the plots to avoid proliferation of figures). 

In general, the results for the NUHM are similar to those of the CMSSM, but preferring 
lower values of mo, which slightly improves the chances of detection at LHC. 

6. Can still SUSY be natural? The impact of the Higgs discovery in the 
naturalness of SUSY 

The results of sets. 4, 5 show that the CMSSM and NUHM scenarios require large super- 
symmetric masses, which challenges the natural coexistence of SUSY with the electroweak 
scale (the problem is actually quite general for the MSSM, as discussed below). There are 
more-or-less sophisticated criteria in the literature to quantify the amount of fine-tuning 
of a given model [89, 90, 91, 92, 93, 94, 95, 96, 97, 98]. However, the most rigorous and 
free-of-prejudice way to treat the fine-tuning is precisely through the Bayesian analysis, 
which automatically incorporates a fine-tuning penalization, encoded in the Jacobian fac- 
tor of eqs.(2.11, 2.13). On the other hand, it is important to note that a Bayesian analysis 
starts with the assumption that the scenario under consideration is the true one, and then 
it gives the regions of higher probability in the associated parameter space. It does not 
directly give a hint about the plausibility of the scenario, unless one establish a comparison 
with an al ternative one. 

In spite of this, it is very interesting the fact that we can use the Bayesian approach 
to study how the plausibility of the MSSM has changed after the Higgs discovery, as we 
are about to see in this section. For that purpose we will use a similar strategy to the one 
followed in subsect.4.3. to study the tension posed by the g — 2 data. As in that case, we 
separate the experimental data in two subsets 



{data} = {@,D} 
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Figure 10: ID marginalized posterior probability distribution of the NUHM parameters: 
mo, ma, mi/2, A and tan/3. The red(cyan) corresponds to flat (logarithmic) priors. 



where T> represents now the Higgs mass data and D the rest of experimental information. 
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Figure 11: 2D marginalized posterior probability distribution of the NUHM in m 1( / 2 -moi m-H — mo 
and rrig — rriq planes. Left panels show the flat priors case and right panels the logarithmic one. 
The inner and outer contours enclose respective 68% and 95% joint regions. 



As explained in subsect. 4.3 we construct now the ££— parameter, defined by 

p{@ ohs \D) p(£i ohs ,D) 
p(@ max \D) ~ p(f max , D) 
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Figure 12: As in fig. 10 but with an additional constraint from DM (WMAP 7-ycars and 
XENON100) considering that the LSP is the soley DM component. Here the cyan line corresponds 
to a positive fj, whereas the blue dashed to a negative //. 



Here f^ max is the would-have-been most probable value of the Higgs mass in the MSSM 
context, given the values of the other experimental data. Our analysis shows that, restrict- 
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Figure 13: As in fig. 11 but with an additional constraint from DM (WMAP 7- years and 
XENON100) considering that the LSP is the soley DM component. 



ing ourselves to the CMSSM scenario, this maximizing value is rrih ~ 117 GeV. In other 
words this would have been the central value of the CMSSM 'prediction" for rrih just before 
the LHC discovery. This can also be seen from ref. [27]. Then the ££— parameter defined 
above gives the relative global likelihood of the actually measured with respect to the 
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maximal one. The results of our computation of Jz? are 

In ££ = —1.92 (log priors) 

In if = -0.70 (flat priors) (6.3) 

which correspond to 1.9a and 1.18cr respectively. (For this computation we have not in- 
cluded dark matter constraints.) The previous results show that, although the experimental 
Higgs mass is not certainly the preferred value in the CMSSM context, it is not as improb- 
able as it would seem. This is in part due to the fact that an important CMSSM region, 
from the point of view of statistical weight, is the Focus-Point one. This region corresponds 
to large scalar masses (but with low fine-tuning due to the small value of fi) and thus to 
naturally large Higgs mass. The results for the NUHM are similar. 

One can summarize this statement by saying that the experimental value of the Higgs 
mass is at ~ 2 a off the supersymmetric predicted central value. Certainly, as past expe- 
riences showed, a 2 a discrepancy is not particularly dramatic and physicists are hardly 
impressed by it. Notice however, that this implies that, if the CMSSM is true, it lies in the 
~ 5% less probable region, which is roughly consistent with the estimates in the literature 
that the 'supersymmetric fine-tuning is of a few percent" . 

So the conclusion is that, although the CMSSM or the NUHM are quite unlikely to 
show up at the LHC (at least using the present strategies to find out SUSY), they are still 
interesting and plausible models after the Higgs discovery. Let us say they are at ~ 2 a. 

Finally, we would like to address a question related to the naturalness of SUSY, namely: 
To which extent the problems of the CMSSM and NUHM remain in general MSSMs? To 
address this question, let us recall that the original motivations for the CMSSM are 1) Min- 
imal Flavour and CP violation, 2) Simplicity, 3) The fact that it arises in some theoretically 
motivated scenarios (like minimal SUGRA or Dilaton-dominated SUSY breaking). From 
these motivations only the first one is robust, but the experimental results do not require 
completely universal soft terms. E.g. the third generation of squarks and sleptons could 
have very different masses. The degeneracy of gaugino masses at Mx is not experimentally 
justified either. Therefore, going beyond the CMSSM (or the NUHM) is very plausible. 
Does this solve the problems of the CMSSM? 

This question can be explored by promoting the CMSSM to more general MSSMs 
compatible with the flavour and CP constraints. Playing with the initial parameters it is 
possible to have a lighter third generation, as invoked in "natural SUSY" models. It also 
allows certain types of supersymmetric spectrum that can evade the detection at LHC. In 
particular if the LSP is heavy or the masses of the supersymmetric states are not far from 
each other ( "compressed spectrum" [99] ) , then the supersymmetric events have small prs 
and are triggered out in the experimental analyses. It is undoubtedly possible to arrange 
the MSSM spectrum in these ways in order to fool the LHC. But, in the absence of a solid 
theoretical reason for them, it sounds a rather artificial possibility. 

In any case, in any MSSM the third generation cannot be too light, since it is needed to 
be rather heavy in order to reproduce (thanks to the radiative corrections) the experimental 
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mass of the Higgs. This requirement is alleviated if Xf is close to its maximizing value, 
but still the required stop masses imply that the electroweak breaking is quite fine-tuned, 
although perhaps not as much as it is often said, as we have discussed in this section. This 
situation could only improve substantially if one goes beyond the MSSM, which are out of 
the scope of this paper. 

7. Comparison with previous literature 

There are several works in the literature that have incorporated the impact of a Higgs 
at ~ 125 GeV and/or the new XENON100 data into Bayesian and frequentist analyses 
[6, 7, 8, 11, 9]. Some of them are previous to the announcement of the Higgs discovery, 
and thus have a slightly different central value and uncertainty of nih from that reported 
in ref.[l, 2]. However, that departure is not dramatic, even though small shifts in the Higgs 
mass have an important effect in the SUSY studies, as is clear e.g. from the discussion 
in sect. 4. On the other hand, there are important differences between the approaches 
followed in the various analyses and the results obtained in them. The aim of this section 
is to show the consistency of all the results, to explain the origin of the discrepancies and 
to argue what are the most reliable results and conclusions. 

Let us start with other Bayesian studies of the CMSSM and/or the NUHM after the 
Higgs discovery, namely those of refs.[7, 8, 11]. One can check that, in spite of being based 
on more or less the same data, the results do not look the same (compared with each other 
and with the present paper). To be more concrete, consider the probability distribution 
function in the m^ 2 — itiq plane in the CMSSM with fi > 0, using log priors and including 
the DM constraints. In ref.[8] they allow for multi-component DM, so the top- left plot of 
their figure 1 must be compared to the top- left plot of our fig. 8. Actually, they are quite 
consistent (the preferred region is in the area of the Focus-Point), although the plot of 
ref.[8] is cut at 8 TeV and does not show most of the allowed region visible in our fig. 8. 
This cut is not an artifact of the representation, but the actual upper limit considered for 
rriQ in ref.[8]. But otherwise, the results are in very good agreement. On the other hand, 
the left plot of figure 2 of ref.[7] and the first plot in the middle row of figure 4 of ref.[ll] 
correspond to the requirement that the DM is supersymmetric. Hence, they should be 
compared with the top-left plot of our fig. 5. Actually, the plots of those two references are 
quite consistent with each other, but they show a preferred region, with quite small soft 
masses, which does not show up at all, neither in our plot nor in the plot of ref.[8] (which 
being less restrictive should display it). What is going on? 

The origin of these discrepancies lies on the different ranges used for the soft parame- 
ters. Not only one can leave aside a most interesting region if the ranges are too restrictive. 
In addition, as discussed after eq.(2.3), the results obtained by marginalizing the non-shown 
parameters may depend on the ranges of integration. The upper-limits of the ranges for 
(m , m 1/2 ) in refs. [8, 7, 11] were (8 TeV, 5 TeV), (4TeV, 2 TeV) and (4 TeV, 4 TeV). This 
explains why the latter two references missed the region where the bulk of the probability 
is (the Focus-Point) and show only the less satisfactory stau-coannihilation and A— funnel 
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Log priors 
CMSSM, |i>0 




m 1/2 (TeV) 

Figure 14: 2D marginalized posterior probability distribution of the CMSSM in the roiw-Wo 
plane, using log priors and DM constraints. The yellow and blue areas enclose respective 68% and 
95% joint regions. The white one corresponds to 99.9% c.l. 



regions. As discussed in subject. 4.2 these two regions are disfavored, as they corresponds 
to rather small soft terms, difficult to reconcile with the Higgs mass (they are still con- 
sistent thanks to the theoretical uncertainty assumed in the supersymmetric calculation). 
This is why these regions do not belong to the 68% and 95% c. 1. regions in the plots of ref. 
[8] and the present paper. In fact, one can force these two regions to appear by plotting 
the contours of higher confidence levels. We have checked that they become visible for the 
99.9% c.l. region, as fig. 14 shows. The new contours are consistent with those shown in 
refs. [7, 11], but they are extremely disfavored. 

The conclusion here is that it is very convenient to use an approach that makes the 
actual ranges for the parameters irrelevant. This is in fact possible, taking profit of the 
Jacobian factor discussed in sect. 2, which penalizes fine-tuned regions at high scale. In 
this sense, we honestly think that the results presented here are very reliable and represent 
an improvement with respect to previous analyses. 

Let us finally comment briefly on frequentist approaches. These are based on the 
analysis of the likelihood function in the parameter space. The scan of the parameter space, 
evaluating the likelihood based on the \ 2 i leads to zones of estimated probability (inside 
contours of constant x 2 ) around the best fit points. In doing so, the unplotted variables 
are optimized to obtain the best x 2 • Two recent frequentist analyses incorporating a Higgs 
mass at 125 GeV are those of ref. [9] and ref. [11] (the first one was made before the Higgs 
discovery). They are quite consistent with each other, obtaining a best-fit point in the 
stau-coannihilation region. The preferred areas are inside this region and in the A— funnel 
one. Roughly speaking, they coincide with those obtained in the Bayesian analyses of 
refs. [7, 11]. The fact that the Focus- Point region with pure Higgsino LSP does not appear 
is probably a consequence of the limited ranges used in the scan of mj^, rriQ. In addition, a 
frequentist analysis is not sensitive to the fine-tuning. Fine-tuning has to do with statistical 
weight and a frequentist analysis is based entirely on likelihood, i.e. the ability to reproduce 
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the experiment, and thus cannot "see" the fine-tuning. E.g. a point (o a region) in the 
parameter space can present an optimal likelihood, but only after an extreme tuning of 
the unplotted parameters. This feature makes this point extremely unplausible (even if 
it is capable to reproduce the experimental data), but this does not affect its privileged 
condition in a frequentist approach. This fact can favor points in the low-energy regions, 
which are suppressed in the Bayesian approach. 

8. Conclusions 

The most obvious conclusion from the results presented in this paper is that if the CMSSM 
or the NUHM are the model chosen by nature, the Bayesian analysis shows that they 
probably live in a region of the parameter space of large supersymmetric masses (squark and 
gluino masses above ~ 2 TeV), which makes very problematic their discovery at LHC. The 
quite large Higgs mass, ~ 126 GeV, is the main responsible for these negative results. 
Actually, given the experimental value of m/j, it would have been extremely improbable 
that the LHC had already detected a signal of new physics corresponding to the CMSSM 
or the NUHM. The Higgs mass is more powerful than the direct detection limits in pushing 
the MSSM parameter space into the large-energy region. Actually, we will be quite lucky 
if one of these models is correct and is still detectable in the next years, at least with the 
current techniques to look for physics beyond the SM at the LHC. We have quantified these 
probabilities with a rigorous and improved Bayesian analysis, whose results are remarkably 
independent of the choice of the priors and of the ranges assumed for the independent 
parameters. 

The previous conclusion is even more dramatic if one includes dark matter constraints 
in the analysis, in particular the last XENON100 data. These are able to discard a large 
portion of the Focus-Point region, while other regions (co-annihilation and A— funnel) are 
very disfavored from the experimental Higgs mass. There survives a fraction of the Focus- 
Point region with large scalar and gaugino masses, where the LSP neutralino is almost 
pure Higgsino. Its mass is constrained to be quite close to 1 TeV. This extreme conclusion 
is relaxed if one allows the dark matter to have other components (e.g. axions). Then a 
large part of the Focus-Point region is rescued. Interestingly, the future XENON1T data 
will be able to test these regions. Therefore, if the CMSSM or the NUHM are the true 
models of nature, it is more plausible to discover them in experiments of direct-detection 
of dark matter than in the LHC. 

Our Bayesian analysis contains an automatic (non-imposed) penalization of fine-tuned 
regions. This penalization arises from the Bayesian analysis itself when the value of Mz 
is treated on the same foot as the rest of the experimental information (and not as a 
constraint of the model). This fine-tuning penalization is encoded in a model-independent 
Jacobian factor, which should never be ignored. This factor is also the responsible for the 
great stability against changes in the ranges of the parameters. We have shown that the 
different choices made for these ranges in previous analyses (which did not include the 
Jacobian penalization) have led to apparent discrepancies between the results obtained by 
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different groups. We demonstrate that all of them are consistent once the freedom in the 
ranges is taken into account. This has motivated also that some previous analyses have 
missed important regions in the supersymmetric parameter space. 

Although the 125 GeV Higgs pushes SUSY to regions where is hardly detectable at the 
LHC, this does not necessarily mean that SUSY is very fine-tuned. Actually, the Focus- 
Point region is not particularly fine-tuned, despite having large supersymmetric masses, 
since the value of \i is small. In this paper we have evaluated to which extent the Higgs 
discovery has diminished the naturalness of SUSY (focussing on the CMSSM for this ex- 
ercise). It is important to note that a Bayesian analysis starts with the assumption that 
the scenario under consideration is the true one, and then it gives the regions of higher 
probability in the associated parameter space. It does not directly give a hint about the 
plausibility of the scenario, unless one establish a comparison with an alternative one. Still, 
it is remarkable that we can use the Bayesian approach to study how the plausibility of the 
MSSM has changed after the Higgs discovery. The key is to perform the same Bayesian 
analysis with the same data except the Higgs mass, say as if we were in the year 2011. 
Then one obtains a preferred value for the Higgs mass (namely, around 117 GeV) with 
a probability profile; and one can compute how far apart from this prediction the actual 
experimental value is. This can be done with the help of Bayesian tecnhiques (comparing 
the evidence for different values of m^). The result is that the Higgs mass is at ~ 2 a 
off the supersymmetric predicted central value. Certainly, as past experiences showed, a 
2 a discrepancy is not particularly dramatic and physicists are hardly impressed by it. 
Notice however, that this implies that, if the CMSSM is true, it lies in the ~ 5% less 
probable region, which is roughly consistent with the estimates in the literature that the 
'supersymmetric fine-tuning is of a few percent" . 

So the final conclusion is that, although the CMSSM or the NUHM are quite unlikely 
to show up at the LHC (at least using the present strategies to find out SUSY), they are still 
interesting and plausible models after the Higgs discovery; and the chances of discovering 
them in future experiments of direct detection of dark matter are quite high. 



A. Focus-Point behavior in the CMSSM and the NUHM 



As first noticed by Feng et al [100], the focus point is a scale at which, for fixed m.^, "io 5 
Aq and tan 8 in the CMSSM, the H u mass, m 2 Hu , does not depend on the initial (high- 
energy) value of rriQ. The value of that scale, which denote Qfp, depends mainly on the 
top yukawa and thus on tan/3. Surprisingly Qfp — Mz- 

Let us now consider the required value of \x to obtain a correct electroweak breaking. 
First notice that from the minimization of the (tree-level) scalar potential, one gets 

„ m 2 iT — m 2 TT tan 2 B \ „ 

Therefore, for large values of tan (3 (> 8) , 

f^-m 2 Hu . (A.2) 
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Consequently if we are in a regime where m\ is small at low-energy and tan /3 > 8, then 
the value of fi will also be small. [The precise value of rn? H ^ depends on the CMSSM 
parameters (mo, mi/2, Aq, tan/3).] All this has important implications for the amount 
of fine-tuning (the smaller fi, the less fine-tuning) and for dark matter (if fi is small the 
lightest neutralino is mostly Higgsino, and the second neutralino and the lightest chargino 
are quasi-degenerate). On the other hand, the scale at which the electroweak breaking is 
evaluated, say QeWj does not coincide exactly with Mz (and thus with Qfp). Rather, it 
is given by an average of the stop masses, which is larger than Mz- The actual value of 
m 2 H u at Qew is given by the running from Qpp to Qew- This represents a positive shift. 
Since m 2 ^ is normally negative at low energy, its absolute value decreases. And from 
eq.(A.2), the value of /x 2 decreases too. This is interesting, but if Qfp is too high (which 
typically occurs for large m 1 / 2 and/or mo), fi 2 becomes negative. This means that there is 
no solution for the electroweak breaking. This puts upper limits on the values of m 1( / 2 and 
niQ (more precisely on a certain combination of them). Near that limit fi 2 is small and we 
can have a pure Higgsino LSP. This is the CMSSM region favored in the Bayesian analysis 
when one includes dark matter constraints. Notice that it takes place for large m-i/2 and 
m . 

Let us now go to the NUHM. In this case, one has the freedom to vary mo and m# at 
high energy in an independent way. One can check that if one leaves rriH(Mx) fixed while 
increases (decreases) mo(Mx), then the value of rn 2 H ^ at low-energy will decrease (increase). 
Consequently, one can decrease the absolute value of mf^ , and thus of \j? , simply by going 
to a region in the parameter space where mo < m#. This region is analogous to the 
Focus- Point region of the CMSSM; but now one can have smaller values of mo, and thus of 
squark masses. Actually, the Bayesian analysis selects this region as the preferred one in 
the NUHM when dark matter constraints are included. This explains the (non-dramatic) 
differences between the CMSSM and the NUHM plots (e.g. Fig. 5 vs Fig. 18). 
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